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Abstract 

Loop calculations involve the evaluation of divergent integrals. Usually one computes 
them in a number of dimensions different than four where the integral is convergent and 
then one performs the analytical continuation and considers the Laurent expansion in powers 
of e = n — 4. In this paper we discuss a method to extract directly all coefficients of this 
expansion by means of concrete and well defined integrals in a five dimensional space. We by- 
pass the formal and symbolic procedure of analytic continuation; instead we can numerically 
compute the integrals to extract directly both the coefficient of the pole 1/e and the finite 
part. 
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Introduction 



Feynman diagrams computations are often a lengthy and hard task. The complexity of typical 
computations grows as a factorial with the number of external legs and/or loops involved in 
the processes. In the Standard Model the main difficulties arise from QCD corrections and the 
production of jets in the final states; but it is not difficult to imagine extensions of the Standard 
Model, in which this problem could seriously limit our capabilities of studying and understanding 
new physical phenomena; in particular if they are involved (in addition to ) by new rather 
strong couplings, which give rise to multiparticle/multiloop amplitudes. This issue motivates a 
rich and increasing activity. Considerable progresses have been achieved in the past 0: new 
methods inspired by string theory [3-5] , helicity amplitudes [6-9] , different recursive algorithms to 
calculate QCD dual amplitudes [|10||, or the scattering matrix elements of generic processes with 
arbitrary initial/final states[|l2]. Very useful approximations have been proposed [12-16], when 



the exact matrix elements are unknown. In loop computations interesting simplifications [17-26] 
can be used. 

While analytical methods can provide us in some specific processes with very simple results, 
a general approach for generic lagrangian and final states cannot avoid the use of the computer 
power. A possible strategy is to implement the Feynman rules into some automatic code with 
the help of some familiar packages for symbolic manipulation of algebraic formulae. However this 
way to proceed becomes ineffective, when the computer has to manage very complex formulae (as 
often it is the case). 

This suggests a second strategy which tries exploiting the computer power with pure numerical 
algorithms and thus avoiding lengthy symbolic manipulation. In the recent past these algorithms 
have been successfully used for the computation of several tree level amplitudes both in QCDl 



2^ and in electroweak processes [29-31]. In this paper we discuss the possibihty to apply pure 



numerical techniques [^] in loop calculations. 

The problem of performing loop computations can be divided in two parts. The first part 



is similar to tree level computations [jTT[: from a given lagrangian one has to produce an explicit 
function of virtual and real momenta which is equivalent to the sum of all the Feynman diagrams 
contributing to the process. Once this task is accomplished, loops calculations have an additional 
problem with respect to tree level ones: integrals over virtual loop momenta are affected by ultra- 
violet and/or infrared divergences. These divergences require the introduction of a regularization 
procedure which usually obliges us to perform all virtual loop integrations in an analytical and 
essentially symbolic fashion (while in tree level calculations one usually performs the integration 
over the real momenta by numerical montecarlo methods). This is a serious difficulty in any 
numerical approach. Here we address this last issue. We will define a numerical procedure to 
perform virtual loop integrations even when these include divergences. 

Let us consider a function f[lfi,Pfi] of the external momenta and of the virtual momenta 
The n dimensional integration [|^] is performed in a region of complex values of n where the integral 
is convergent; then one considers the analytical continuation to values of n where the integral is 
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not well defined. In general for n ~ 4 the analytical continuation can have a pole Xje (e = n — 4). 
We can use the Laurent expansion around e ~ to define some functionals lk\j\ such thatQ 

I d"//[/^] = + /o[/] + h\]V + h\fY + • • • (1) 

Since the n dimensional integration is well defined for any function /, also the functionals lk\!\ 
are unambiguously defined from equation (|1|). For instance, /o[/] coincides with the ordinary four 
dimensional integral 

/ d'm,\ (2) 

if / is a convergent function. While the left-hand side of the (0) is essentially a formal definition 
which does not lead us to an obvious numerical integration procedure, each Ik\f\ can be written 
in terms of well defined and convergent integrals. In the following we show how to build explicit 
integral representations for the functionals In the next sections, we explain two different 

definitions (even if equivalent) of the /a,.[/]; in the final part some very simple numerical examples 
will be given in order to emphasize the practical purposes of the method. 



Method of Integration A 

Let us assume that we have to perform an integral in n > 4 dimensions. The virtual momentum 
has n components. We call the components with /x < 4 and the components with /x > 4. 
Then we can rewrite the integral 

I ^^if%,v,\ = I d'mi^p,] dn = I d'if[i,,p,f~' d/ da (3) 

d^e is the solid angle in the subspace (/i > 4) of dimension e orthogonal to the four- vector l^. 
Note that is invariant under rotations in this subspace, since the external momenta do 

not have components /x > 4. We can exploit this invariance as follows. First we rotate the vector 
L such that h ^ Q and Z„ = for any n > 5. Thus we can write In = (liyhJsyh, Vi) with 

2 2 2 

t = ih) + ih) + ■■■ + (In) and omit all the components with /i > 5. We have reduced our n 
dimensional space into a five dimensional space. Second, we can perform the integral in dfle and 
we obtain from eq. (H) 

5f(/d*VP„pj)*"-dt (4) 
Now the integral in dt can be integrated by parts 

d'l^,f[l„P,]) dt. (5) 



ei2 + e)\J dt 



^At one loop. 
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The function / depends on t through the fifth component of l^. The function t^+i has a cut in 
the positive real axis of t ( Re{t) > and Im{t) = ) , in fact the rotation t — >■ te*^"^ gives 
— ). ;fi+2+*'r^. Therefore we can apply the following identity 



oo 1 

where the contour of the last integration in t must contain all the poles and singularities lying 
in the complex t plane and must avoid the cut on the positive real axis (see figure 1). Then the 
integral (H) above becomes 

1 7r^/2 r 4t^+t / r \ 

^r3i=)Fi/i(2T7)(/'"i?«"'"'0^'- « 

Now we can make the expansion in powers of e (the Euler constant and log(7r) have been re- 



absorbed into a redefinition of the scale /i as in the MS) and we obtain the final result 

/^(^ + (-i + iogM])) (/ d'l^fii^p,]) dt. 



In the above expression each integral is now convergent, due to the action of the derivative d^/ dt^ 
appearing inside the integral in d^l. Thus one can separately calculate the singular and the non 
singular part, simply evaluating well defined and convergent integrals. We can check the result 
above in a specific example and see how the above formula can be used for practical calculations. 

Suppose that we have to evaluate the following divergent integral 

d'^^TT^- (9) 



P + m 



In our approach is not a n dimensional object but it is a more concrete five dimensional vector 
with components (Ji,l2,h,h, V^); thus P = l^l^ = P + t with P equal to the squared length of 
the real four- vector {h,l2,hyU) and t a complex number. If we are interested in the finite part, 
we can apply the non singular term in equation (Rl) and we have 



it , . , , / d^ P + t 



-(-i + iogHD^^P^^j d'/d. (10) 

After the derivative the integral in d'^l is convergentf], and this yields 

/|(-i + iogNl)(-^)d* (u) 

The integral over the contour in the complex plane is equal to the residue at the pole t ~ — 

-7r^(- l + log[m2])m^ (12) 



^We consider only value of t with non zero imaginary part. See the discussion at the end of this section. 
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Our final result is correct. The advantage of the procedure above is clear: eqs.(p!0|-p!l|) are well 
defined and concrete expressions, any step can be done in a pure numerical way. Note that 
the integral in d^I must be done after the derivative d^/ dt^, and before the integral in dt; 
otherwise the integral in d^I would be non convergent. In practice, if the above procedure is done 
numerically, one should replace the integral in dt with a sum over a finite set of points tf, for the 
convergence of all integrals, it is enough to choose the ti in this set with a non zero imaginary 
part. 



Method of Integration B 

The formula (|^) is a transparent and compact expression, which makes manifest some properties 
of the dimensional regularization; for instance, the invariance under the translations ^ 
is obvious, since it comes directly from the translational invariance of the integral in d^I, which 
is convergent (after the derivative in d^/ dt^). However in certain numerical calculations the 
(H) could be not efficient. We discuss a different method which is more powerful in practical 
calculations, even if it requires a slightly more involved formula. 

There are several different ways of rewriting the (^); each of them depending on the way we 
choose to parameterize the virtual momentum Here we do not intend to make an exhaustive 
study, we will simply describe a quite general procedure to obtain well defined integral represen- 
tations for the Ik in (|l|). 

First we observe that each Ik[f] is a linear operator 

4[/i + /2] =4[/i] + 4[/2] (13) 

4 [A/] = Xhif]. (14) 

It is quite natural to think that the linearity (|13|-0) implies that the 4 are authentic integrals or 
sum of integrals. In fact we can introduce a simple trick to build concrete integral representations 
for generic linear functionals. 

Consider the space of functions which admits a Laurent expansion around a point to- They 
are defined by a set of through the expansion]^ 

oo 

m= E «n(t-to)". (15) 

n=— oo 

A linear functional /[/] acting on this set of functions can be written 



/[/] 



J2 a„(t-to)" = E a„/[(t - to)"] = E "-nbn (16) 



^We also assume that the series converges strongly enough to justify the steps below. 
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where = I[{t — to)"] is a real (or complex) number. Then we can use the identity 

1 {t-t,T 



C 2TTi {t - to) 



k+1 



dt = 6nk (17) 



for any n,k = 0, ±1, ±2, ... if the integral contour C is a small circle, in the complex plane of t, 
around the point to- Using the (|1^, the (|T6D can be written 

h, 

dt = 



-| „ oo oo 

/[/] = ^JiT. -nit-tor){ E 



^ f fmt] dt 



k=^oo (t — to) 



k+1 



with 



This is an integral representation of the linear functional /. This simple trick can be used to build 
explicit and compact integral representations for any linear functional, whose action on each term 
of a Laurent expansion is known. 

Let us apply it to our problem. We start defining the angular integration in n = A + e 
dimensions. The analytical integration of any tensor with k (even) indices and constructed with 
the n dimensional vector yields^ 

J IfilpU ■ ■ ■ Ipl-ph dl^4+e = ^^^ (2 + ! + £-)!! ''^^^^^^ • • • + • • • + 9!iag-pD ■ ■ ■ g-pa) ^^^''^^^ (20) 

dn^+e = sin2+^[7] . . . sin2[^]sin[0] d7 d^ d0 d5 (0 < 7, ^, . . . , < vr; < 5 < 27r) is the solid 
angle. It is understood that the indices {ft, u . . .) above are contracted with some external momenta 
and/or other tensorial indices (like the spin of the external particles), since the scattering matrix 
is lorentz invariant. Thus0 only indices with /i < 4 are relevant in our problem. This has been 
emphasized by the small bar on the Greek indices. 

We look for an integral definition which reproduces the (|20|). To begin the procedure we 
parameterize 1^ with an array of five components^ 

III ^ t {x cos[6'], X sin[6] cos[0], a; sin[6'] sin[0] cos [5], x sin [6'] sin[0] sin [5], Vl — x'^^ (21) 

9 (f) and 5 are the three angles in the space of four dimensions, t and x are two complex variables 
which will be integrated over a complex contour as explained below. The need of the auxiliary 

'*Here we assume that the overah factor 7r5/r[-|] has been re-absorbed mto a redefinition of the scale /i, as 
prescribed by the MS scheme. 

^We remind that only has components > 4, and factors l^l^ = are not relevant for the angular integration. 

^In practice this means that any vector (including external momenta etc.), must be written as an array of five 
components, but we keep in mind that only has a non zero fifth component. 
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variable x will be clear in a moment. Suppose that we integrate over the four dimensional solid 
angle dfl^ = sin^[6']sin[0] d9 dcj) d6 all the tensors built with the array (^TD in place of an n 
dimensional Z^. We get 

r (2)" 

J Iplpla . . . Iplpla dVLi = 271^ ^ 'j^y^^ igppgpa ...go,p + ... + g^agpp ■ ■ ■ gpa)t^x^- (22) 



Comparing this result with the ( poD we see that tensorial structure is identical, but the overall 
factor is not correct. The term of order e° can be obtained simply setting x = 1 and = t^. 
Instead for the term of order e, it would be enough to replace x^ with a suitable factor^ 



-^h = ^-{l-^E-i'[2 + k/2]) (23) 



for any even k. 



This can be achieved by means of a linear functional I such that I[x^] = bk- In fact, following 
the recipe above (eqs.(|T5|)-(p!9D) we can build the integral below 



nf] = — <fy:. -fe. m dx = — / + ~ /m dx (24) 

where the integral contour is a circle around the singularity x ~ 0. The integrand above has a 
cut in the segment of the real axis between -1 and 1. It is easy to see that in the limit of a path 
of integration very close to this cut, the logarithm can be approximated by its imaginary part. 
The real part does not contribute to the full integral, and the logarithm can be replaced by izvr. 
Taking also into account that f[x] = f[—x] we can rewrite the ( p^ 

nf]=e ['(t^) f[x]dx. (25) 



V 1 — 



Here the notation (...), means that 



' ih[x])J[x] dx = r h[x] if[x] - /[I]) dx, (26) 
10 Jo 

and this makes the integral convergent near x ~ 1. 

Finally we are able to write a compact formula for the angular integration in 4 + e dimensions 

J /[g dn,+, = ( J /[/,] dn,)^^^ + KiV ^^^'^ (t^)+ "^"^ '^^') + ' ^^^^ 

The last step is to define the integration in dP = dt^. Again we follow the trick (|T5| ) - ([T9|) : 
we have to compute the b„ in order to get the function w[t'^]. By definition, we know that for any 
a>2 



^ The function V' is the weh known derivative d/ dx log[r[x]]; it comes out when we expand the factorial in 
the (pQ) and take only the first order in e. 
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(-2 + a)(-l + (T) (m2) 

for a=l,2 



^ + ^(---), (28) 
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61 = 2 h m^log[m^] + £(. . .) 

62 = =-^-(log[m2] + !)+£(...) (29) 

and for cr=0,-l,-2,... the integral is zero {b„ = 0). If the integration contour in the t'^ complex 
plane is a small circle around the pole singularity = — and replacing the 6„ in the (|19D with 
the 60- above we get 



[f] = E b.+i{t' + m^)'^ = -f(^- + log[- t'] ) (30) 



00 

r ,9n 

(T=0 



and 



/ dtV[t1 = -^(^ / dt^) - / log [ - ^1 dt^) + • •)• (31) 

The integral contour is closed, it must contain all the singularities in the complex plane of the 
function /[t^], except for the singularity in = and the cut for real and positive in all 
integrals where the logarithmic function appears (see for example figure 1). We will comment on 
this contour of integration later on. 

Combing the angular integration ( P7| ) and the ( pij ) we get the full formula 

J /[/,] d"/ = J dn4+. = ^-f + Io + he + ... (32) 



where 

/^i[/[y = (i//[^/^]^'dtMfi4)^^^ (33) 

(34) 

Clearly the integrals above are well defined and can be evaluated numerically.^ For example to 
evaluate the first integral of Iq in the expression (^4]), the contour of integration in dt^ can be 
chosen as in figure 1. Namely we can divide the contour in three paths. The path A is very close 

^ Here all integrals are understood to be in Euclidean space. In some kinematical regions, the integration (even 
if convergent) needs to be regularized by the prescription — > + ie. Also for some infrared singularities the 



34|) must be rearranged differently or one should use the method A. The practical use of these methods in various 
realistic situations certainly demands a much more extensive discussion, in this letter we simply present the general 
idea, a more complete and systematic study will be done elsewhere (see ref . ) . 
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to the cut of the logarithm. This function can be replaced with iivr (the sign depends if we are 
above or below the cut). This yields the following integral 

Pi = dt' dn,) . (35) 

^ •'^ir x = l 

Note that this is an ordinary four dimensional integral with an infrared and an ultraviolet cut-off: 



Im r * 




Re t' 



Figure 1: Integration contour in the complex plane. The thick line in the positive real axis 
represents the cut of the logarithmic function. The path A is very close to this cut and the 
logarithm can be replaced there by its imaginary part. and Af^ are respectively the infrared 
and ultraviolet cut-off. 



in fact after setting a; = 1, the array (^) may be regarded as a genuine four- vector. Then we have 
the paths B and C, that are two circles with radius and A]^ respectively 



^ £ f[i,]t' log [ - 1^] dt^ dn4)^_^ + £ f[i,]t' log [ - 1^] dt^ da 



x=l 



(36) 



These integrals can be seen as some counterterms which cancel the Auv and Air dependence of the 
integral (|35|). The sum is equivalent to the finite part ( in the limit e— »-0) of the dimensionally 
regularized integral. 

Through the linearity (|T3|-|I^), this method can be generalized to two (or more) loop calcula- 
tions. If and are the two loop virtual momenta, then we can rewrite the integration variables 
(Tl = /2{i+e/2) (iri4+£/2 and d"g = dg,, g^^^+'^Z^) dq^ dfis+j. g,, is the component of g^ parallel 
to g^ is the total length of the components of g^ orthogonal to Z^. One then defines two six 
dimensional vectors and g^, analogously to (pi|), as functions of five angles 9,(f),6,(l)',S', two 
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auxiliary variables x and y (needed for the integrations dfl^^^ and dl^a+g) and the complex vari- 
ables Z^, Qj^. Applying the trick (|15D-(|1^) one obtains the function w[l^, q^^, q^]. Needless to 
say, this contains dilogarithms in addition to some logarithms. From this guidelines, one gets the 



analogue of the (34), for two loop calculations [33 



Some Numerical Examples 

In order to make clear the practical purposes of these methods, we discuss here some very simple 
examples. Suppose that we want to integrate 

(Tl \ (37) 

{I + p) + ra^ 

with rn? = 2 and = 1. Our numerical approach must reproduce the analytical result 



2 



2— + ttW \og[m^] + ... (38) 



e 



Suppose that we want to extract the finite partQ vr^m^ log[m^]. We apply the second method B 
described above. We use the five dimensional array ( pTl) for 1^. Then the integrand is a function of 
X, t^, 9, (p, 6. The finite part of the ( pTj) is obtained computing the two integrals in Jq (eq.(0)). For 
the first one we choose the contour of integration as in figure 1: the numerical values of A^^, and Air 
have to be chosen in such a way that the contour A+B+C contains all the physical singularities 
in the complex plane. Since the integral is infrared convergent we can take the limit Aj^ — > 0. 
Instead, in the ultraviolet region, we can cut the integral at A^^ = 9. We also must set x = 1. 
Then the integral (p5|) becomes 

= r dt^ dn^-- ^— = 51.929. (39) 



2t2 + 2tcos[^] +3 

One can recognize in the integral above the ordinary four dimensional integration with a cut-off 
A^„. Then the integration in dt^ follows the path B 

P2 = — / dt'^ dni\og\~t^]-- ^— = -43.1815. (40) 

2m Jb ^ ^2t2 + 2tcos[^] +3 ^ ^ 

This completes the first integral in Iq. The second one in eq.(|^) contains a non trivial integration 
in the variable x, which must be performed using the prescription ( [26| ) . The contour of integration 
in dt'^ is rather simple: there is no logarithm, and no cut in the real axis; then we do not need 
to follow the path A. There is no infrared singularity and also the path C vanishes. The only non 
trivial contribution comes from the path B, a circle of radius A^^. This gives 

P^ = — t i ^—Tw, de( dx dn^ = 4.935. (41) 

i27rio + 2xtcos[g] + 3 \l - x^ J + ^ ' 



^ It is more instructive to discuss the finite part, since the integral contour in (^) is trivial, and extracting the 
singular part is straightforward. 
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Note in the denominator the appearance of the variable x, which comes out when we take the 
scalar product / ■ p with from eq.(pl|) and p'^ in the direction fi = 1. The sum of the three 
contributions yields 

^o[- k ;] =Pi + P2 + P3 = 13.6825 ~ vrWlog^], (42) 

[l + p) + 

in perfect agreement with the analytical expression. The method can also be applied for infrared 
divergent integrals. The integral 

r 1 1 

/ 2 (43) 

is ultraviolet convergent and infrared divergent. Thus we choose A^^ = oo and Aj,. = 1/4. We get 

roo 1 11 

P^= e dt^- dn- = 10.843, (44) 

7l/4 2 t2 + 2tcos[^] +3t6 ' ^ ' 

Po = — i dt^ dnAos\-t^]^- ^— = -12.1255 



2m Jc 4 OL J f2 + 2t cos[e] + 3 
Pi = — { dt^f dx d^]4 = -0.1825. (45) 

The sum of the three integrals gives 

^°[^77— 4— ^] = Pi + P2 + P3 = -1.465, (46) 
i {l-\-p) 



very close to the exact result 

- 27- V* ' ^ 4 



-i-vr^fl + logft). (47) 



Conclusions 



In the dimensional regularization of loop integrals, usually one defines the integral as a function 
of the number of dimensions in a region of n where the integral is convergent. Then one makes 
an analytic continuation to obtain a definition of the integral also in regions of n where the 
integral is divergent. This procedure is clear and unambiguous, however it can only be applied 
in pure analytic (and symbolic) calculations. It cannot be converted into an obvious numerical 
procedure. One is obliged to perform symbolic calculations, which sometimes becomes lengthy 
and/or unaffordable. In this paper we have shown how to set up a different approach, where each 
coefficient the Laurent expansion in powers of e in eq. (|I]) can be written in terms of concrete and 
convergent integrals (well) defined in a five dimensional space (see eq.(||) and eqs.(^),(^3D and 
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(|3^ ) . Instead of abstract n dimensional vectors we have to deal with "concrete" five dimensional 
vectors, and we have to perform integrations over a compact space. This formulation has the 
advantage to allow us the evaluation of all integrals through pure numerical methods. Clearly 
further work is needed, to generalize the above result to more loops and to prove the efficiency 
of this technique in realistic physical problems. However we believe that the simplicity of the 
numerical approach makes it very promising. 
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